integrating factor method for pdes